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Abstract 

Using the machinery of smooth scaling and coarse-graining of observ- 
ables, developed recently in the context of so-called fluctuation operators 
(originally developed by Verbeure et al), we extend this approach to a rig- 
orous renormalisation group analysis of the critical regime. The approach 
is completely general, encompassing classical, quantum, discrete and con- 
tinuous systems. Our central theme is the analysis of the famous 'scaling 
hypothesis', that is, we make a general investigation under what cluster 
conditions of the 1-point correlation functions a scale invariant (non-trivial) 
limit theory can be actually attained. 



1 Introduction 



One of the crucial ideas of the renormalization group analysis of, say, the critical 
regime, is scale invariance of the system in the scaling limit. This is the famous 
scaling hypothesis (as to the underlying working philosophy compare any good 
text book of the subject matter like e.g. Q and references therein). Central in 
this approach is the socalled blockspin transformation j^j. That is, observables 
are averaged and appropriately renormalized over blocks of increasing size. At 
each intermediate scale a new effective theory is constructed and the art consists 
of choosing (or rather: calculating) the critical scaling exponents, so that the 
sequence of effective theories converge to a (scale invariant) limit theory, provided 
that the start theory lay on the critical submanifold in the (in general infinite 
dimensional) parameter space of theories or Hamiltonians. 

Usually the calculations can only be performed in an approximative way, the 
main tools being of a perturbative character and are frequently model dependent. 
Typically, in the more general discussion, spin systems are employed, to motivate 
and explain the calculational steps. While the general working philosophy, based 
on the concepts of asymptotic scale invariance, correlation length and the like, 
is the result of a deep physical analysis of the phenomena, there is, on the other 
side, no abundance of both rigorous and model independent results. 

This applies in particular to the control of the convergence of the scaled l- 
point correlation functions to their respective limits if we start from a microscopic 
theory, lying on the critical submanifold. In this case, correlations are typically 
long-ranged and the usual heuristic arguments about the interplay between poor 
clustering, on the one side, and formation of block variables of increasing size, 
on the other side, become rather obscure as one is usually cavalier as to the 
interchange of various limit procedures. It is known, that this is a very dangerous 
attitude in this regime. 

Furthermore, the clustering of the higher correlation functions in the various 
channels of phase space may be quite complex and non-uniform in general. A 
concise and selfcontained discussion of the more general aspects and problems, 
lurking in the background together with a series of note and references, can be 
found in ||, section 7. 

Usually, the crucial scaling relation (the scaling hypothesis) 

W?{Lx x ,... ,Lx f ,ix*) =L- ln -L l ^-W?{x 1 ,... ,x l]f i*) (1) 

which is conjectured to hold at the fixed point (denoted by \x* in the parameter 
space), is the starting point (or physical input) of the analysis. Here, Wf denote 
the truncated Z-point functions (see below), L is the diameter of the blocks, Lxi 
are the respective centers of the blocks, n is the space dimension, 7 the statistical 
renormalisation exponent. If it is different from n/2, we have an 'anomalous' 
scale dimension. 
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In the following analysis we want to concentrate on the derivation of such (and 
similar) scaling relations for the /-point functions from the underlying microscopic 
characteristics of the theory. We will do this in a completely general way, that is, 
the underlying model theory can be classical or quantum, discrete or continuous. 
We try to adopt a quite general point of view, making as few model assumptions 
as possible. Our strategy is it, to deal only with the really characteristic (model 
independent) aspects of the subject matter. 

In this sense our approach is perhaps similar to the one, expounded in e.g. [|J 
in the context of the analysis of the ultraviolet behavior in algebraic quantum field 
theory. A nice review of the unifying aspects of the cluster of ideas,underlying 
renormalisation with applications to numerous fields is also ||. 

At the end of this introduction we want to remark that there exists a superfi- 
cially different approach, which is more related to the well established concepts of 
renormalisation theory in quantum field theory (see, for example, |J, || or [J7J), 
that is, renormalising propagators, Greens's functions and path integrals and in 
which scale invariance is present on a more implicit level. 

We do not intend (and actually feel unable) to relate all these different aspects 
to each other in this paper. We will instead briefly sketch, what we are going 
to do in the follow sections. We concentrate our analysis entirely on the hierar- 
chy of correlation functions which define the theory. We construct renormalized 
limit correlation functions from them which happen to be scale invariant, thus 
defining a new limit theory. These limits are set into relation to the degree of 
clustering of the original microscopic correlation functions. We do not openly 
discuss the flow of, say, the renormalized Hamiltonians through parameter space. 
The characteristics of these intermediate theories are however implicitly given by 
their hierarchy of correlation functions as was explained in e.g. p) or H. 

The same holds for the effective Hamiltonians which emerge on the intermedi- 
ate scales in the ordinary approach. In our approach the effective time evolution 
is carried over from the microscopic theory as described in H or (in a slightly 



other context) in [[|, see also [IT]. In case we work in an environment, defined by 



ordinary Gibbs states, our framework would exactly yield these effective Hamil- 
tonians. In our paper we treated the infrared limit. As our scaling approach 
depends on a continuous parameter, R, which could as well go to zero, a slight 
extension of our approach would allow to study also the ultra violet behavior. 
In this case we have to work, however, in a space-time picture, that is, the time 
coordinates have also to be scaled. The cluster properties are then expected to 
be even more intricate as the clustering with respect to both space and time is 
of a markedly dynamical nature. 
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2 Concepts and Tools 



As to the general framework we refer the reader to |§ . One of our technical tools 
is a modified (smoothed) version of averaging (such a possibility is also briefly 
mentioned in the notes in 0). Instead of averaging over blocks with a sharp cut 
off, we employ a smoothed averaging with smooth, positive functions of the type 





for \x\ 


1° 


for \x\ 



f R (x):=f(\x\/R) with f(s) = I £ ~ (2) 



(one can choose even more general classes of functions). This class of scaled func- 
tions has a much nicer behavior under Fourier transformation, as, for example, 
functions with a sharp cut off, the main reason being that the tails are now also 
scaled. We have 

f R (k) = const ■ R n ■ f(R ■ k) (3) 

Remark: One might perhaps think that this choice of averaging will lead to a 
different limit theory. This is however not the case. Furthermore, the mathe- 
matical differences between the two approaches, that is, using sharp or smooth 
and scaled cut off functions, are relatively subtle and not so apparent. We will 
investigate these aspects in |10| . 



Another point, worth mentioning, are the implications of translation imparl- 
ance. We have for the correlation functions 



W(xt, ... , x n ) = W(xi - x 2 , ■ ■ ■ , x n -i - x n ) (4) 

The truncated correlation functions are defined inductively as follows (see ||) 

w( xu , x n ) = n ■ ■ ■ > x ^ ( 5 ) 

part P i 

The (distributional) Fourier transform reads 

W T (p 1 , . . . ,pi) = W T (p 1 ,p 1 +p 2 , . . ■ ,Pi H Pi-i) ■ 5{p x H pi) (6) 

The dual sets of variables are 

i 

Hi := Xi - x i+ i , qi = ^2pj i<(l-l) (7) 
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3 The case of Normal Fluctuations 

As in 0, we assume that away from the critical point the truncated /-point 
functions are integrable, i.e. G L 1 (i?™^ -1 ^), in the difference variables, yi := 
Xj — Xi + \. As observables we choose the translates 

A R ( ai ), A R ( ai ) , A R (a) := R~ n/2 ■ J A(x + a)f(x/R)d n x (8) 

(where, for convenience, the labels 1 . . . / denote also possibly different observ- 
ables). We then get (for the calculational details see |§, the hat denotes Fourier 
transform, translation invariance is assumed throughout, the const may change 
during the calculation but contains only uninteresting numerical factors): 

(Ar(oi) • • • A R (ai)) T = const ■ R ln / 2 ■ 
J f(Rpi) ■ ■ ■K-R\px+- ■ ■+Pi-i})-W T (p 1 , . . . ^.e-^-WgfciErp* J[d Pi 

= const ■ R ln/2 ■ R-V- 1 ^ ■ 

J M) • ••/(-&+• • ■+p1 1 })-w t ( p ' 1 /r, . . . ,p{_ 1 /i2). C -Si-V^-.. e *-«E!rV < /H JJ dp j 

(9) 

We now scale the a^'s like 

Oj :=R-Xi , Xi fixed (10) 

This yields 
{A R {R-X l )---A R {R-X l )) T = 

const ■ R?- l >' 2 ■ J e-^i-V,* • e lX ^^ ■ 

/Vi) ■ ■■f(-[Pi + • • • + p'i-i\) ■ W T (p[/R, . . . ,pli/R) U d Pi ( n ) 

As the / are of strong decrease and W T continuous and bounded by assump- 
tion, we can perform the limit R — > oo under the integral and get: 

Case 1 (I > 3): 

Km (A F R (R ■ X,) ■ ■ ■ A F R (R ■ X t )) T = (12) 

K— *oo 

Case 2 (/ = 2): 

lim (^(i? • • X 2 )) T = canst ■ [ W T (0) ■ e ~*i ■ f(p[) • f(-p'iWi 

R^oo J 

(13) 
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Conclusion 3.1 In the normal regime, away from the critical point, where we 
assumed L 1 -clustering, all the truncated correlation functions vanish in the limit 
R — > oo apart from the 2-point function. We hence have a quasi free theory in 
the limit as described in or in the work of Verbeure et al ( cf. the references in 



4 The Relation to the Heuristic Scaling 
Hypothesis 

In the following sections we will develop two rigorous, technically slightly differ- 
ent, approaches which implement the physically well- motivated but, nevertheless, 
to some extent heuristic scaling hypothesis. The analysis will be performed in 
coordinate space, Fourier space, respectively. In this section we restrict ourselves 
to the two-point correlation function, for which the asymptotic behavior is con- 
siderably simpler and more transparent. 

The pedagogical reason for representing two technical methods is the follow- 
ing. It turns out that, as may be expected on physical grounds, critical systems 
and their characteristics are sitting exactly on points in the parameter space 
where also the necessary mathematical methods turn out to be of a singular 
character. This can best be exhibited by contrasting the two methods developed 
below. 

Remark: In the rest of the paper we replace the exponent n/2 in the definition 
of Ar(o) by an exponent 7, which will usually be fixed during or at the end of a 
calculation. It plays the role of a critical scaling exponent. 

Let us hence study the behavior of 

(A R (R ■ X X )A R {R ■ X 2 )) T = R~*r ■ J W T (( Xl - x 2 ) + R{X X - X 2 )) 

■ f(xi/R)f (x 2 / R) dxidx 2 
= iT 27+2n J W T {R[{xi - x 2 ) + (X x - X 2 )\) ■ f(x 1 )f{x 2 )dx 1 dx 2 (14) 

We make the physically well-motivated assumption that-in the critical regime- 
W T decays asymptotically like some inverse power, i.e. 

W T (xi -x 2 ) ~ F(x! -x 2 ) ■ \xx -x 2 r {n ~ Q) 0<a<n (15) 

for \x\ — x 2 \ — ► 00, F bounded and well-behaved. 

From the last line of ([14]) we see that, as / has compact support, we can 
replace W T , for (Xi — X 2 ) 7^ and R — > 00 by its asymptotic expression and 
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absorb the remaining (xi, x 2 ) -integration in an appropriate constant. We hence 
get for R large: 

{A F R (R ■ X 1 )A^(R ■ X 2 )) T w const ■ R- 2 ^ 2n . R-("-«) . \x x - X 2 \^ n ~ a) (16) 
Choosing now 

7 = (n + a)/2 (17) 
we get a limiting behavior (for R — > 00) as 

~ \X ± -X 2 \- {n - a) (18) 

Central in the renormalisation group idea is that systems on the critical sur- 
face (i.e., critical systems) are driven towards a fixed point, representing a com- 
pletely scale invariant theory. This idea is usually formulated in an abstract 
parameter space of, say, Hamiltonians. In our correlation function approach the 
fixed point shows its existence via the scaling properties of the correlation func- 
tions, that is 

(L • (X - Y); //) = L- 2(n ~ 7) Wf (X - Y; //) (19) 

with fi* describing the fixed point in the (usually) infinite dimensional parameter 
space. We see from the above that this is exactly implemented by our limiting 
correlation functions, as we have (with the choice 7 = (n + ot)/2): 

W^X - Y; a*) ~\X- Y\- [n - a) (20) 

That is, the above scaling limit leads to a limit (i.e. fixed point) theory, exactly 
reproducing the asymptotic behavior of the original (microscopic) theory. But 
there is yet another important point worth to be mentioned. 

In the more general situation of /-point correlation functions we do not neces- 
sarily have such a simple decay behavior. Furthermore, it would be desirable to 
have precise information on thresholds up to which a certain method leads to rig- 
orous and reliable results. We therefore describe and illustrate two methods with 
the help of the relatively transparent behavior of the 2-point functions, which 
will then be applied to the more complex and general situation of the higher cor- 
relation functions. Using the assumed simple asymptotic behavior of the 2-point 
functions, we proceed now as follows. 

4.1 Method One 

We assume the existence of a certain exponent, a, so that (x 2 denoting the 
vector-norm squared) 

G{x) := W T (x) ■ (1 + x 2 Y n ~ a)/2 = const + F(x) (21) 



7 



with a decaying (non-singular) F. Fourier transformation then yields: 



R~ 2J - J W?(( Xl - x 2 ) + R(X 1 - X 2 ))f{x 1 /R)f{x 2 /R)dx 1 dx 2 

= R^- J G((x 1 -x 2 ) + R(X 1 -X 2 ))>[l + ((x 1 -x 2 ) + R(X 1 -X 2 ))^™V 2 - 

f(x 1 /R)f(x 2 /R)dx 1 dx 2 



/ 



e - iRp{xi - X2) ( R - 2 + _ x j + (Xi _ x 2 )) 2 )-("-«)/ 2 /(a;i)/(x2)dxida; 



(22) 



where we made the substitution x — > R ■ x. 

We now assume the support of / to be contained in a sufficiently small ball 
around zero (or, alternatively, (Xi — X 2 ) sufficiently large so that (xi — x 2 ) + 
(X 1 — X 2 ) 7^ for Xi in the support of /). With 

G(p) = const ■ S(p) + F(p) (23) 

the leading part in the scaling limit R — > oo is the 5-term. Asymptotically we 
hence get for R — > oo (setting y := x\ — x 2 Y := X-y — X 2 ): 

Rn +a~2 1 . const • y |y + y | "("-«) • / * /(y)dy (24) 

with 

f*f(y):= J f(y + x 2 )-f(x 2 )dx 2 (25) 

The reason why the contribution, coming from F(p), can be neglected for R — > 
oo is the following: / is assumed to be in X>; by assumption the prefactor never 
vanishes on the support of f(xi). Hence the whole integrand in the expression in 
square brackets is again in V and therefore its Fourier transform is in S, that is, 
of rapid decrease. We can therefore perform the i?-limit under the integral and 
get a rapid vanishing of the corresponding contribution in R for each p ^ 0. As 
the 5-contribution has been extracted from G(p), this proves the statement. 

As / * / has again a compact support, we have that, choosing 

7 = (n + a)/2 (26) 

the limit correlation function behaves as ~ \Xi — X 2 |~( ra ~ a ) as in the above 
heuristic analysis. 



8 



4.2 Method Two 

As in the case of normal clustering or last section), one can, on the other 
hand, improve the too weak decay of W T {x\ — x-z) and transform it into an 
integrable (i.e. L 1 -) function. So, with a similar notation as in the preceding 
subsection, we choose a suitable exponent a in 

P a {x l -x 2 ):={l + \x l -x 2 \ 2 ) a ' 2 (27) 

so that 

G(y) := W T (y) ■ P' 1 G L 1 (y := Xl - x 2 ) (28) 

In contrast to Method One, there is of course a whole range of such possible 
exponents, a > so that 

G<y) = ll» T a> """ < 29 > 

I f. L tor a < ctinf 
Proceding as in Method One, we get 

iT 27 J W T {y + R-Y)- f(x 1 /R)f(x 2 /R)dx 1 dx 2 = 

R . 2l _ Rn+a J dp Q( p / R ^ . e -i P Y . J e -wy( R ~2 + ( y + Y ) 2 ) a l 2 ■ f * f(y)dy 

(30) 

Again the obvious strategy seems to be to choose 

1 = {n + a)/2 (31) 

and perform the limit R — ► oo. With the same support properties as above, that 
is, y + Y 7^ for a;i, x 2 G support of /, the integrand in square brackets is again 
infinitely differentiable with respect to y. Hence, its Fourier transform is fast 
decaying in p. 

Remark: Note that for a/2 non-integer and without the above support restriction, 
there will show up a singularity in sufficiently high orders of differentiation for 
vanishing R~ 2 . One can however control these singularities and show that the 
analysis still goes through in the case where the support condition does not hold. 
One gets however some mild constraint on the admissible a's. 

Therefore we can again apply Lebesgues' theorem of dominated convergence and 
perform the i?-limit under the integral. This yields the expression 

G(0) • j dpe' ipY ■ J e -*» .\ y + Y\ a -f* f(y)dy = 

const ■ G(0) • J ' 8{y + Y) -\y + Y\ a ■ f * f{y)dy = const ■ G(0) • (32) 
result of the above support condition). 
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Conclusion 4.1 With a chosen so that G(y) G L 1 and 7 = (n + a)/2, the 
limit can be carried out under the integral and yields the result zero. This shows 
a fortiori that there is no a m i n with the property that there is a non-vanishing 
limit-two-point function. Put differently, we have an cti n f but no a m i n (cf. 

So, in contrast to Method One, the relevant exponent, a in f, is of such a pe- 
culiar nature that we definitely cannot apply the above method of interchange of 
taking the limit R — > 00 and integration. But nevertheless, we will show that 

7 := (n + a inf )/2 (33) 

is the correct critical scaling exponent leading to a sensible limit theory and that 
this ai n f is exactly the a, we have determined in Method One. 

We have learned above that in order to arrive at a non-zero limit correlation 
function, we are definitely forbidden to exploit Lebesgues' theorem of dominated 
convergence in the above expression. The reason for the vanishing of the respec- 
tive expression was that with 

lim G{p/R) = 6(0) (34) 

R^oo 

we have to evaluate f g(p)dp with 

g(p) := J e-^y +Y ^\y + Y\ a -f* f(y)dy (35) 

This integral happens to be zero due to the explicit factor, \y + Y\ a and the 
assumed support properties. 

So, we have to investigate what happens for a = a in f. As we learned above 
that there is no a min , we can conclude 

Observation 4.2 For a = cti n f, G a (y) is no longer in L l , with 

G a (y):=W T (y)-(l + y 2 r a / 2 (36) 

We know that for a < cti n f the decay of G a (y) is so weak that the Fourier 
transform develops a power law singularity in p = 0; that is, we can conclude 

Lemma 4.3 For ai n f — a := e, G a has a singularity of the form \p\~ e nearp = 0. 
For a = ai n f the singularity is of logarithmic type near p = 0. 

This statement can again be proved by a scaling argument. Let G a have a non- 
integrable tail of the form \y\~^ n ~ £ \ For the (distributional) Fourier transform 
we then have 

G a (\ ■ p) = const ■ J e iXpy ■ G a {y)dy = const ■ X' n ■ J e ipy ' ■ G a (y'/\)dy' (37) 
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For A — > we can, as above, replace G a by its asymptotic expression, which goes 
as \y\~^ n ~ £ ^ and conclude that G a (Xp) contains a leading singular contribution 
~ A _e (modulo logarithmic terms). We hence see that 

G a (p) ~ (38) 

near p = as a distribution (that is, the above resoning is to be understood 
modulo the smearing with appropriate test functions; see e.g. [0). For a = otj n /, 
the singularity must be weaker than any power, that is, must be of logarithmic 
type. 

By Method One we get a limit correlation function which clusters as \Y |-( n_a ). 
One may wonder where this behavior is hidden if we use Method Two. Taking 
only the singular term in G(p/R) into account, we have (with 7 := (a in f + n)/2) 

lim R^{A R {RX 1 ) ■ B R (RX 2 )) T ~ lim const ■ [ \n(\p\/R) ■ g{p)dp (39) 



and g(p) as in equation (35). We can again neglect the term 

InR- g(p)dp (40) 



as J g(p)dp = 0. 

Assuming at the moment that a were an integer (we will get the general 
result by a scaling argument), the prefactor \y + Y\ a can be transformed into 
corresponding ^-differentiations of f * f(p), which, by partial integration, can 
then be shifted to corresponding differentiations of ln(|p|). This transformation 
yields an expression of the type \p\~ a times a smooth and decaying function. 
That means, we essentially end up with an expression like 



dp \p\ a ■ e 



-ipY 



e~ tpy f * f(y)dy 



(41) 



By the same reasoning as above we conclude that the singularity, \p\~ a , goes over, 
via Fourier transform, into a weak decay proportional to \X% — X 2 \~ <yn ~ a \ that is, 
we arrive at the same result as in Method One, whereas the reasoning is a little 
bit more tricky. 

For a general non-integer a the argument could be made precise by analysing 
the distributional character of an expression like r 13 with r := \x\ and its Fourier 
transform. As the analysis is a little bit tedious, we refer the reader to [ |l^] . 
On the other hand, one can use a scaling argument as above (with Y := A ■ 
Yq , Yq fixed as A — > 00). This yields an asymptotic behavior of the form 



A~ (n " Q) • J dp ln(|p|) • J e ~ lXp{y+Yo) ■ \y + Y Q \ a ■ f * f(Xy)dy 



(42) 
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The evaluation of the integral for A — > oo can be done as follows: As / * / 
has compact support, the volume of the support of f*f(Xy) shrinks proportional 
to A~ n . Therefore the expression in square brackets scales as ~ A _ra . On the 
other hand (due to an 'uncertainty principle' argument), its essential p-support 
increases proportional to A n . That is, the two effects compensate each other and 
we have again a large- F behavior ~ |y|-( n_Q ) as before. 

We conclude that both methods lead to the same aymptotic scaling behavior 
of the renormalized two-point function. 



5 The General Cluster- Analysis at the Critical 
Point 

We now study the general situation of the presence of some long-range correlations 
in the 1-point functions. In contrast to the much simpler situation prevailing 
in the case of two-point functions, the clustering may be quite complicated, in 
particular, the dependence on the number, I, i.e. the number of observables, 
occurring in the expressions, may be non-trivial. Therefore, we have to investigate 
these aspects in more detail. 

From general principles (see e.g. [13|]) we know that in a pure phase there is 
always a certain degree of clustering. We make the slightly stronger assumption 
that it is in some way of the kind of an inverse power law at infinity (to be 
specified below). We want to study the scaling limit of 

(A R {R-X 1 )---A R {R-X l )) T (43) 



with 



A R (a) : = iT 7 J A(x + a)f{x/R)d n x (44) 



and an, at the moment, unspecified exponent, 7. 
The above expression can be written as 



» 1 1 

/ W T (( Xl - x 2 ) + R{X 1 - X 2 ), . . . , - x t ) + R(X^ - X,)) • Y[f(xi/R) ]Jdx 

J i=l t=l 

(45) 

Fourier transformation yields (with W T (qi, . . . , qi-i) considered as a distribution 
on S(R( l -V n )) 

const ■ R 1 ^ ■ j W T ( gi , ... , ft _ a ) ■ e-^ R ^ ■ 



i=l i=l 



/L 1 L b 

e-i^RPM . e iR*-i*i Y[f(xi) Y[dx 



i=i i=i 



1-1 



Y[d Pj (46) 
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with Yj := Xj —Xj + i and the wellknown relation between the q- variables and the 
p- variables (see e.g. section 2 or ||). For calculational or notational convenience 
we will employ both sets of variables which are linear combinations of each other. 

As / is in T> by assumption, the Fourier transform of Yl f( x i) is m $ an d the 
function in square brackets is a function of (Rpi, . . . , Rpi-i) or (Rqi, . . . , Rqi-i), 
being of rapid decrease in either set of variables. As a consequence, for R — > oo 
and at least one pj being different from zero, the expression approaches zero faster 
than any inverse power (together with all its derivatives). 

From this we see that, as R — > oo, the region of possible singular behavior 
is located around (p)!" 1 = or (q) 1 ^ 1 = 0, implying also pi = — X)i 1 Pj = 0- 
We can hence infer that only the singular behavior of W T in (q) = will matter 
in this limit. As a consequence, it will be our strategy to isolate this singular 
contribution in W T and transform it in a certain explicit scaling behavior in R, 
which can be encoded in some power, R~ a , in front of the integral. 

The singular behavior of W T (q) at (q) = is related to the weak decay of 
W T (y) at infinity. The limiting behavior of W T (y) can, however, not expected 
to be simple or uniform (at least not in the generic case) as (yi, ■ ■ ■ ,yi-i) or 
(xi, ... ,xi) can move to infinity in many different ways. We may, for example, 
have that (xi) together with all \x{ — Xj\ go to infinity or, on the other side, the 
variables move to infinity in certain fixed clusters of finite diameter. The rate 
of decay of W T {y) should of course depend in general on these details. Corre- 
spondingly, the singular behavior of W T (q) in the infinitesimal neighborhood of 
(q) = should depend on the direction in which (q) = is approached, that is, 
the limit may be direction-dependent. 

In the light of this general situation we must at first decide, in which kind 
of limit we are mainly interested. Inspecting the expression (f|3|), we actually 
started from, we choose in a first step our fixed vectors, (Aj), so that 

Xi- Xj^ for alii, j (47) 

As a consequence, all distances, \RXi — RXj\, go to infinity for R — ► oo. As 
in the preceding section, we can choose the support of / so small that, with 
G supp(f), we have 

\R(Xi-Xi) -R{X j -x j )\ -> oo (48) 

In this particular case we may expect a relatively uniform limit behavior on 
physical grounds. 

Remark: Similar problems occur in quantum mechanical scattering theory. 

Under this proviso the following assumption seems to be reasonable. 

Assumption 5.1 Under the assumption, being made above, we assume the fol- 
lowing decomposition ofW^(y) to be valid: It exists a function, (1 + H(y)), H{y) 
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homogeneous and positive for y ^ so that 

G(y) : = (1 + H(y)) ■ W T (y) = const + F(y) (49) 

with F sufficiently decaying at infinity in the channel, indicated above, i.e. {\y%\ — > 
oo for alii = 1, . . . ,1 — 1} and 



H{Ry) = R a 'i ■ H{y) 

Remark 5.2 A typical example for H{y) is yf ) a ^ 2 ■ 
Fourier transforming G(y), we get 

G(q) = const ■ 6(q) + F(q) 



(50) 



(51) 



and expression ([I6| ) becomes (compare the related expression in Method One of 
the preceding section) 



const ■ R l{n ~^ ■ ! Yl dp 3 G{q) ■ e"^^ • 
l 

/i i 
e -i£ l f 1 RpjXj . e i^i-m . ^ + H ( Ry + BY))- 1 . Yl fix,) ■ JJ dxj 
i i 



(52) 



By assumption, H is homogeneous of degree a[. So we can extract a negative 
power of R, R~ ai , from the expression in square brackets. Furthermore, we 
observed above that for R — > oo only the vicinity of q = matters. Finally, by 
assumption, the contribution coming from F(q) can be neglected in this limit 



(compare the corresponding discussion in the subsection [4~T| ; as a consequence 
of the assumed support properties, the expression in square brackets is again 
strongly decreasing). We hence have 

. i-i 

lim (A R (R ■ X x ) ■ ■ ■ A R (R ■ Xi)) T = lim const ■ #(^-<*Ht) . / TT d Qj 

J 1 

^ i=l i=l 

-iZi^RpjXj . e iR*-ixi . (R-<*'i + H(y + Y)y ■ Y[f(xi)Y[dxi 

i=i i=i 



S(q)-e 



Remark: We see again the reason for the special choice being made above as to 
the support properties of the functions f(xi), leading to the result yj + Y 3 ■ ^ 
on the support of /. Without this assumption, we see for our above example, 
H{y) = (S yf) a ^' 2 5 that m the limit, where R~ a 'i vanishes, we would get a singular 
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contribution at points where y + Y = in the integrand in square brackets. These 
terms would make the following discussion much more tedious. 

If we now make the choice 

7 := 7i = n- a'Jl (54) 
we arrive at a finite limit expression, depending on the coordinates (Xj): 

/i i 
H^y + Y)- 1 ■T\f(x t )T\dx t 
i i 

(55) 

which is a function of the coarse grained difference coordinates 

Yj = X 3 - X j+1 (56) 

For the Yj sufficiently large, it is approximately a function 

W limit (Y) « const -H^Y)- 1 (57) 

That is, the renormalized limit correlation functions reproduce the asymptotic 
power law behavior of the original microscopic correlation functions. 
For later use we introduce the new scaling exponent, ai, via 

a[ + on = (I - l)n (58) 

This implies 

7, = (n + a,)/* (59) 
The underlying reason for this choice is that an asymptotic decay, ~ r -(i-i)n^ j g 

1 /2 

just the threshold for Wj being integrable or non-integrable (with r := yj) ). 

We have arrived at the following result: We are interested in a scaling-limit 
theory for R — > oo. In order to get a non- vanishing and finite limit theory, we 
have to choose the scaling exponent for / = 2 as 

7 = 7 2 = (n + a 2 )/2 (60) 

Furthermore, we have extracted the asymptotic form from the higher truncated 
/-point functions, Wf(y), and have absorbed it in an explicit scaling factor, R 
to some power. If the limit theory is to be finite, the corresponding scaling 
exponents for I > 2 have to be less or equal to zero. This yields unique 7;'s as 
threshold values. 

A corner stone of the philosophy of the renormalisation group is that the 
scaling exponents of the scaled observables remain the same, irrespectively of the 
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degree of the correlation functions in which they occur. That is, these exponents 
are fixed by the exponent, 72, and we have 

7 = 72 > li (61) 

(the latter exponent being derived from equation (j59|)), in order that the limit 
correlation functions remain finite. 

Conclusion 5.3 We have the following alternatives for R —>■ 00: 

72 > 7/ W* R (62) 

72 — 7/ =^ h m is finite and non-trivial (63) 

72 < 7/ => W^ R -> 00 (64) 

72 > 7/ /or a// Z > 3 ; t/ie fixed point is gaussian or trivial. The limit theory is 
quasi-free. The limit theory is non-trivial if ^2 = li for at least some I > 3. For 
72 < 7/ for some I, the limit theory does not exist. 



Remark 5.4 The corresponding analysis can also be done by employing Method 
Two (discussed in the preceding section). One can even omit the support con- 
ditions assumed above. The treatment then becomes more involved but the end 
result is the same. We discuss one particular case below. 

To complete the scaling and/or cluster analysis of the truncated correlation 
functions, we have to analyze the other channels and the respective consequences 
for scaling exponents and cluster assumptions. 

We mentioned several times that without the support condition 

{Xi - Xj) + (xi - Xj ) + (65) 

for Xij G supp(f), the analysis would become more tedious. On the other side, 
this assumption is violated if the observables move to spatial infinity in certain 
clusters. The extreme case occurs when all Xi are chosen to be zero, i.e: 

(A R (l)---A R (l)) T , R^oo (66) 

(the indices 1, ... ,1 denote the different observables). This scenario was already 
briefly discussed in section 7 of in connection with phase transitions and/or 
spontaneous symmetry breaking, which are also typically related to poor spatial 
clustering. 

With the same notations as above we have 
(A R (1) ...A R (l)) T = const -R 1 ^ ■ 



i=l i=l 



/t, — u c — V 

e-iEi-'Jfcyay . e ^-i** . JJ f( Xt )Y[dx 
i=i i=i 



1-1 



Y[d Pj (67) 
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Assuming again the existence of a suitable homogeneous function, Hi(y), in this 
channel, we get asymptotically two contributions 



/i i 
Hi{y)- 1 -\{f{x l )\{dx % 
i i 



(68) 



and 



i-i 



canst ■ R l ^)-«'i • j J] d Qj F(q u qi _ x ) 



i=l 



■1 i 



\\dxi 



i=i 



(69) 



The first term has almost the same form as above. But now the function in 
square brackets in the second contribution is no longer of strong decrease as the 
integrand (considered as a function of (x)or (y)) is no longer in T> as it will have 
a singularity in y = 0. We can however provide the following estimate on the 
degree of this singularity of Hf 1 in y = 0. We assumed throughout in this section 
that Wj is not integrable at infinity, that is, the clustering is weak. On the other 
side, this asymptotic behavior is exactly encoded in H^ 1 , as we observed above. 
The threshold where integrability goes over into non-integrability for H^ 1 is a 
behavior 



' i-i 



We can therefore conclude that 

a[<(l- l)n 



(70) 



(71) 



in the above construction if Wf is non-integrable at infinity. If a\ is even strictly 
smaller than (/ — l)n, which is the ordinary case in the critical region, we have 



Observation 5.5 

a[ < (I - l)n 

implies that H^ 1 is integrable near y = 0. Hence 

i 

1 

due to the compact support of f . 



(72) 



(73) 
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From this we infer again that, with 



7; = n - oc'Jl = (n + a L )/l (74) 

the contribution (^) is finite in the scaling limit. For the contribution ( |6"§D we 
have by the same reasoning that the function in square brackets is a continuous 
function of (Rq), which goes to zero for Rq — > oo (due to the Riemann-Lebesgue 
lemma) . 

On the other side, we have no precise apriori information about F(y) and 
F(q). F(y) goes to zero at infinity as the asymptotic behavior is contained in 
if -1 , but its rate of vanishing is not clear. 

Conclusion 5.6 If the integrand of contribution (j^) is lying in some LP , so 
that the limit, R —>■ oo, can be performed under the integral, the whole expression 
vanishes in the scaling limit. 

In this situation we are left with again with the first term, which is the limit of 

„ i i 

const • / Hfa + F) 1 J] f(x t ) J] dx t (75) 
J i i 

for Y — > 0. That is, in this case it holds 

Theorem 5.7 If the situation is as in the conclusion, Wl im (Xi, . . . ,Xi) is con- 
tinuous and we have in particular 

m (0, . . . , 0) = lim W} im (X h ... } X t ) (76) 

We can hence resume our findings as follows: If the assumptions, made above, 
are fulfilled and if the functions, Hi, can be chosen consistently in all channels, 
so that the 7z's, resulting from the relation 

7z = (n + an) /I (77) 

are smaller than or identical to 72, we arrive at a full limit theory, being well- 
defined in all channels. In this case the renormalization group program works 
and yields a non-trivial scaling limit. 
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